Does a Surface Polariton Have Spin? 
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We consider a p-polarized surface electromagnetic wave (a classical surface polariton) at the 
interface between the vacuum and a metal or left-handed medium. We show that the evanescent 
electromagnetic waves forming the surface polariton inevitably possess a backward spin energy flow, 
which, together with a superluminal orbital energy flow, form the total Poynting vector. This 
spin energy flow generates a well-defined (but not quantized) spin angular momentum of surface 
polaritons which is orthogonal to the propagation direction. The spin of evanescent waves arises from 
the imaginary longitudinal component of the electric field which makes the polarization effectively 
elliptical in the propagation plane. We also examine the connection between the spin and chirality 
of evanescent modes. 
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Introduction. — The spin angular momentum (AM) of 
light arises from the circular polarization of propagating 
electromagnetic waves and is directed along the wave mo- 
mentum [l[ . It obviously vanishes for a linearly-polarized 
wave. Furthermore, the spin AM is known to be purely 
intrinsic, i.e., independent of the coordinate origin [2|- 
5]. At the same time, propagating light can also posses 
orbital AM: it originates from phase gradients and can 
have both intrinsic and extrinsic parts [H-Q. The spin 
and orbital AM are produced, respectively, by the lo- 
cal spin and orbital energy flows (EFs) [fHa] . These are 
associated with vector and scalar properties of the field 
and together constitute the Poynting vector (momentum 
density). The separation of the spin and orbital parts of 
the AM and Poynting vector is unique, and both parts 
are separately observable - e.g., via the motion of test 
particles 0, d, 1(J- HI or the evolution of instantaneous 
distributions of the wave field [l3| . The properties of the 
spin and orbital EFs in optical fields were recently exam- 
ined in detail apart from those in evanescent waves. 

In this paper we investigate the spin and orbital 
properties of ^mearZy-polarized evanescent electromag- 
netic waves by considering a p-polarized surface polari- 
ton at the interface between the vacuum and a negative- 
permittivity medium [l4- 17|. We demonstrate that, 
despite its linear polarization, the evanescent wave in- 
evitably carries non-zero spin EF and spin AM, the lat- 
ter being directed orthogonally to the wave momentum. 
Moreover, the orbital EF is sitper luminal, whereas the 
spin EF is backward, which together ensures sitfeluminal 
local energy transport in the forward direction. The spin 
of the evanescent wave arises from the imaginary longi- 
tudinal electric field, which generates a rotation of the 
electric-field vector within the propagation plane. Fur- 
thermore j_we examine the relations between spin and 
chirality [HI, [l9j] for evanescent waves. 

Scalar and vector features of evanescent waves. — We 
consider a p-polarized surface polariton plane wave at 



the z = interface between the vacuum (2 > 0) and a 
medium (z < 0) with real permittivity e = e m < and 
permeability fi — fi m . Assuming that the surface mode 
propagates along the x-axis, its unit-amplitude electric 
and magnetic complex fields can be written as T3- Toj] 
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where the "+" and "- 
ties in the z > and 



superscripts denote quanti- 
< half-spaces, and = 



exp [ik p x =p k z — ioJot] are the scalar wave functions lo- 
calized at the interface. Here loq is the frequency, fco = 
loq/c, whereas the evanescent waves / are characterized 
by complex wave vectors = fc„x ± inr^i, which satisfy 



the dispersion relations k = k, 



„± 2 _ 



e^ikQ. Us- 



ing the proper boundary conditions at the interface, this 
yields the surface-polariton parameters 14j-4l6jj : 
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We would like to emphasize two important features 
of the solutions (1) and (2). First, the surface polari- 
ton propagates along the x-axis with the wave number 
k p > fco i an( l its phase velocity is v p h = cko/k p < c. At 
the same time, the local energy-transport velocity (which 
in free space becomes the group velocity) can be deter- 
mined using the relativistic relation between the energy 
W and momentum p: p — v g W/c 2 . For the scalar evanes- 
cent waves f^, the x-component of the momentum and 
energy are proportional to k p and u>o, and we arrive at the 
superluminal group velocity v g o — ck p /ko > c. (In con- 
trast, for a plane wave in free space, propagating at some 
angle with respect to the x-axis, with k x < fco, we would 
obtain w p h = cko/k x > c and v g o = ck x /ko < c.) Thus, 
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FIG. 1: (color online). Instantaneous distributions of the 
real electric field £ (r, t) and intensity T (r, t) for a surface 
polariton (1) propagating along the surface of a metal (z < 0) 
with Em = —1.5 and p m = 1. The electric field in each point 
rotates anticlockwise (clockwise) at z > (z < 0), whereas 
the intensity wave crests move with the phase velocity u p h = 
cko/kp < c. 



The orbital EF is essentially determined by the phase 
gradient of the scalar wave function /, whereas the spin 
EF is produced by the gradients of the polarization ellip- 
ticities cp E = lm (E* x E) and ip H = Im (H* x H). 

The separation (3)-(5) works well in a homogeneous 
medium, but in the presence of inhomogeneities (e.g., 
interfaces), the spin and orbital EFs acquire non-zero 
divergences: V • Ps = — V • Po 7^ 0, which does not 
makes a physical sense. Since V • P = 0, one can modify 
the separation of the spin and orbital EFs, P = P' s + P' , 
such that V • P' s = V • P' Q — (cf. Q). In this manner, 



we obtain P' 



s 



A, P' 



A, with 
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4fc 



[V/U 1 x tp E + Ve 1 x tp H ] 



(6) 



it might seem that the apparent superluminal group ve- 
locity of the scalar evanescent waves contradicts rel- 
ativity. 

Second, consider the polarization of the surface po- 
lariton (1). Although it can be regarded as a linearly - 
polarized p mode with the electric field lying in the prop- 
agation (x, z) plane, we emphasize the imaginary char- 
acter of the longitudinal x-component of the field. It 
arises from the transversality condition E ± ■ = with 
imaginary fc^ = ±ik^. This results in the = F 7r /2 phase 
difference between the E^r and Ef- field components, i.e., 
in the rotation of the electric field in the (x, z) plane. In 
other words, a p-polarized evanescent wave is, in fact, 
elliptically polarized in the propagation plane. 

Figure 1 shows the temporal evolution of the real elec- 
tric field £ (r, t) = RcE (r, t) and instantaneous intensity 
T(r,t) — |ReE(r, t)\ + |ReH(r, t)\ for the surface po- 
lariton (1). The motion of the wave crests demonstrates 
subluminal phase velocity v p h, whereas the electric-field 
vector rotates in each point anticlockwise (clockwise) at 
z > (z < 0). A nice interplay of these features is re- 
vealed below. 

Spin and orbital energy flows. — The time-averaged en- 
ergy density and local EF (the Poynting vector) of an 
electromagnetic wave in an isotropic medium with real e 
and pi are given by [20j 
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cjRe[E*xH], (3) 



where e = d{uje) /duj > 0, jl = d(ujfi) /du > 0, and 
g = (Stt)^ 1 in Gaussian units. The EF determines the 
density of the kinetic momentum of the field, p = P/c 2 . 
Generalizing previous free-space results 0-0] > the Poynt- 
ing vector in the medium can be divided into its spin and 
orbital parts, P = P5 + Po, as 



eg 
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Im [fi- 1 V x (E* x E) + e- 1 V x (H* x H)] , (4) 



eg 

2k 



Im [^- 1 E*-(V)E + £- 1 H*-(V)H] . (5) 



This term describes a "spin-orbit interaction" which van- 
ishes in a homogeneous medium, but becomes important 
at interfaces. 

Due to the above-mentioned polarization properties of 
the p-polarized evanescent waves, the electric-field ellip- 
ticity does not vanish for the surface polariton (1) and 
yields 

<^ = -2fe-- + *y, cp- = . (7) 

Owing to the strong z-gradient, this ellipticity results in 
a non-zero spin EF (4). Substituting Eqs. (1), (2) and 
(7) into Eqs. (4) and (5), we obtain 



p + _ cgn -2k+z- p - 



fcgfcp 



e 2 ^ 2 x, (8) 
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P+ = ^e-^x, P G = -fV 6 " (9) 

Accordin gly, the total Poynting vector of the surface po- 
lariton is |14-[l6t 

P+ = ££*V*« + ** , P- = ^e 2 «^x. 

kp € m kp 



(10) 



Importantly, because of the discontinuity of [i~ 1 ipE at 
the vacuum- medium interface, Eq. (7), strong counter- 
propagating boundary spin and orbital EFs arise there. 
Taking into account the "spin-orbit" correction (6), these 
boundary EFs are 



SP S = SPn = 
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1 - 



1 



S(z). (11) 



Thus, evanescent waves possess a backward spin EF 
(8), Ps || —x, (in the medium it can be forward if jjL m < 
0). This spin EF is subtracted from the forward orbital 
EF (9) to give the total energy current (10). Figure 2 
shows the distributions of the time-averaged field inten- 
sity and EFs (8)— (11) for polaritons on the surface of (a) a 
metal and (b) a "perfect" left-handed medium with e m = 
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Distributions of the time-averaged 
2 , spin (Ps), orbital (Po), and total 



FIG. 2: (color online), 
intensity I = |E| 2 + |H 
(P) EFs (8)-(ll) for the surface polariton (1) propagating 
along the surface of (a) a metal (e m = —1.5, /i m = 1) and 
(b) a "perfect" left-handed medium (e m = /i m = —1). The 
boundary flows (11) are indicated by dashed arrows. In the 
vacuum (z > 0), the backward spin EF is subtracted from 
the forward superluminal orbital EF to provide subluminal 
energy transport. In the metal, the spin EF dominates over 
the orbital EF, which results in the backward energy transport 

El. 



bital AM consists of extrinsic and intrinsic parts 0, 
(L) cxt = (r) x (p ) , (L) hlt = (L)-(L) cxt , 



(13) 



where (r) is the centroid of the beam. 

To prove the intrinsic nature of the spin AM of sur- 
face polaritons, we calculate the integral spin EF (8) and 
(11). Remarkably, the positive boundary flow (11) pre- 
cisely balances the negative bulk flow (8) and (Po) = 
J Ps dz = 0, akin to the case of propagating waves 
Thus, although the spin EF is crucial for the local en- 
ergy transport, it does not transfer energy globally. This 
ensures that (S) cxt = (r) x (p s ) = @H. At the same 
time, the global energy transport is realized by the or- 
bital EF: (P ) = (P) = (cgk /2K+k p ) [l - e m 2 ] x, and 
the ratio (P) / (W) yields the known group velocity of 
the surface polariton [3, [lj| . 

The value of the AM is typically normalized by the 
integral energy (W) [l|. Since (W) is strongly dependent 
on the dispersion in the medium, we first calculate the 
spin and orbital AM for the free-space evanescent field in 
the z > half-space. Using Eqs. (l)-(3) and (8)-(10), its 
energy is (W + ) = J Wdz = g/2n + , whereas the spin 

2>0 

and orbital AM (12) become 



(14) 

Noteworthily, the same spin AM (in units of h per par- 
ticle) can be obtained via calculating the normalized ex- 
pectation value of the quantum spin-1 operator S with 
the fields (1) [8]: 



fi m = — 1 [2l[ . To understand the importance of the spin 
EF, note that Eqs. (8)-(10) in the vacuum can be written 

as P+ = -c (K +2 /k kp) W + ±, Pj = c(k p /k ) TU+x, 

P+ = c (ko/kp) VF+x, where W + = ge~ 2K+z is the en- 
ergy density (3). Using the relation P — pc 2 — v e W, one 
can see that the orbital EF corresponds to the superlumi- 
nal group velocity, v g o = c k p /ko > c, mentioned above, 
while the backward spin EF reduces the total momentum 
and the corresponding group velocity becomes sublumi- 
nal: v g — c ko/kp = Vph < c. Hence, it is the backward 
spin EF that ensures proper local energy transport in 
evanescent electromagnetic fields. 

Angular momenta and chirality. — The spin and orbital 
parts of the wave momentum p = P/c 2 determine the 
spin and orbital AM of the electromagnetic field 
Their spatial densities are given by S = r x p$ and L = 
r x po, whereas the integral (in our case - integrated 
over z) values can be written as 

(S) = (rx P5 ), <L} = (rx Po ). (12) 

The spin represents a purely intrinsic AM, while the or- 



(E+,H+|S|E+,H+) 
(E+,H+ |E+,H+) 

In the whole space, the AM yield 



(15) 



(S) = 1 - 



(S+) , (L) = 1 - 



(16) 

Thus, evanescent waves and surface polaritons posses 
well-defined (but not quantized) spin and orbital AM di- 
rected orthogonally to the propagation (x, z) plane. 

The separation of the intrinsic and extrinsic parts of 
the orbital AM is determined by the centroid of the field, 
(z) = (z W) I (W), which depends on the medium disper- 
sion. As an example we consider polaritons on the sur- 
face of a "perfect" left-handed material with e m (ojq) — 



(wq) 



-1 



2l|. In this case, k + 



the boundary 



EFs (11) vanish, and the EFs arc mirror anti- symmetric 
with respect to the z = plane, Fig. 2b. Choosing the 
model dispersions e m (oj) = \x m (cj) = 1 — 2loq/lo, we have 
£m (k>o) = Am (^o) = 1, and the energy densities become 
mirror symmetric with respect to the z = plane. In 
this case (z) = 0, and purely intrinsic spin and orbital 
AM (16) yield (S) = 2 (S+) and (L) = 2 (L+). 
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FIG. 3: (color online). The same intensity and EF dis- 
tributions as in Fig. 2, but for a surface-polariton wave 
packet propagating along the surface of a dispersive left- 
handed material with electric and magnetic plasma frequen- 
cies u) e — a v /3/2o;o and w M = 2u>o, which correspond to 
£m (wo) = —0.5 and fi m (ujq) = —3. The counter-propagating 
EFs in the vacuum (z > 0) and the medium (z < 0) bring 
about the vortex EFs in the localized wave-packet solutions. 
These solutions carry spin, intrinsic orbital, and extrinsic or- 
bital AM along the j/-axis. 



An intrinsic AM can be associated with a circulating 
EF, i.e., a vortex [US, One can see that the energy 
circulation in the surface-polariton field is non-zero for 
any contour encircling the origin (see Fig. 2). In such 
circumstances, the vortex EF appears upon the longitu- 
dinal ir- localization of the field, i.e., considering a surface- 
polariton wave packet [22j. As an example, we consider 
a realistic left-handed material with plasma dispersions 
e m (uj) = 1 — uj 2 /uj 2 , and fj, m (ui) = 1 — uj 2 /lu 2 . For this 
case, Figure 3 shows a surface-polariton wave packet cal- 
culated numerically using a narrow Gaussian spectrum of 
solutions (1) centered around ui = ujq. It is clearly seen 
that the spin and orbital EFs form counter-circulating 
vortices in the (x, z) plane. In this generic case, the or- 
bital AM (12) and (13) contains both intrinsic and ex- 
trinsic parts. 

Finally, we explore an important connection between 
the spin and chirality of the wave 18 , 19| . Akin to the en- 
ergy density W and energy flow P, one can characterize 
the chirality of the electromagnetic field by the chirality 
density K and chirality flow 4? which satisfy the continu- 
ity equation [3, [3] ■ Generalizing the earlier free-space 
results to the case of a homogeneous medium, these time- 
averaged quantities can be written as 



K = g Im(H*-E) , * = yIm(/i _1 E* x E + e X H* x H) . 

(17) 

Substituting here the surface-polariton field (1), we im- 
mediately arrive at 



K = 0, # = |^V^0. 



(18) 



ellipticity of the field polarization, Eq. (7). To under- 
stand these results, note that for the propagating fields, 
the integral chirality (K) is intimately related to the av- 
eraged helicity of photons, whereas the chirality momen- 
tum (4>/c 2 ) is proportional to the spin AM (S) In 
our case, the helicity vanishes identically because the spin 
AM is orthogonal to the momentum, and (K) = 0. At 
the same time, calculating the integral chirality momen- 
tum, we find that it is indeed proportional to the spin 
AM: 



(*/c 2 ) = 2k (S) . 



(19) 



Thus, the chirality density vanishes since H* • E = 0, 
whereas a non-zero chirality flow is determined by the 



Here we obtained an additional factor of 2 as compared 
to the general result for propagating fields Note- 
worthily, the connection between the chirality and he- 
licity (rather than spin) is quite fundamental. The main 
point is that chirality is a parity-odd but time-even prop- 
erty [23| | . Spin changes its sign upon time inversion, while 
helicity does not. 

Discussion. — We have shown that evanescent waves 
in free space and surface polariton modes at the inter- 
face with a negative-permittivity medium possess super- 
luminal orbital energy flow and non-zero backward spin 
energy flow. The latter originates from the rotation of 
the electric field in the plane of propagation, and it is 
necessary for proper energy transport. The EFs gener- 
ate well-defined spin and orbital angular momenta which 
are orthogonal to the propagation direction of the wave. 
The helicity and chirality density naturally vanish in such 
case. It is worth noticing that the previously consid- 
ered AM of propagating waves [l| and near-field vortices 
24] essentially require the superposition of multiple plane 
waves which produce the necessary gradients. In sharp 
contrast, the transverse spin and orbital AM already ap- 
pear here for a single surface-polariton plane wave (two 
evanescent waves) owing to its natural confinement (in- 
homogeneity) across the interface. 

Our results appeal to experimental tests revealing the 
unusual transverse spin of surface polariton evanescent 
waves. Typically, spin AM manifests itself in interactions 
with probe particles, and it is important to discuss the 
fundamentals of such experiments. 

The spin AM is usually observed in propagating opti- 
cal fields via the spinning motion of the absorbing or bire- 
fringent test particles [J, [lfj| - [l2l | . Any local perturbation 
of the field with a non-zero ellipticity (e.g., a small region 
of the field exclusion around the particle) immediately in- 
duces radial intensity gradients and circulating spin EF 
Q . This circulating EF spins the particle in any point of 
the elliptically-polarized field. In the case of evanescent 
waves, the situation becomes more complicated, because 
any local perturbation in the (2, z) plane will drift along 
the x-axis with velocity v g = v p h = ck$/k p . Still, the 
probing particle can experience a non-zero circulation of 
the spin EF which will induce its spinning motion. 
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In addition to the spinning motion, test particles can 
move linearly in the background EFs. Such interaction 
crucially depends on the physicalproperties of the parti- 
cle. For instance, Berry noticed [8J that the forces acting 
on small absorbing and conducting particles are propor- 
tional to the orbital EF (8) and the total Poynting vector 
(10), respectively. Hence, such particles in the evanescent 
field will experience forces k p /ko times higher and ko/k p 
times weaker than the analogous force from a propagat- 
ing plane wave without spin EF. Thus, monitoring the 
linear motion of different particles, one could observe the 
action of different EFs. 

Finally, the vanishing of the chirality density implies 
that the interaction of the surface-polariton plane waves 
with small chiral particles (e.g., molecules) cannot distin- 
guish between right- and left-handed enantiomers. The 
verification of this conclusion could also be an important 
confirmation of the above theory. 
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